Confrontation of QMD model with the experimental data for 

40^^ +45Sc, i^Mm +197 Au and ^'^^Xe reactions 

Tajinder Pal Singh'* and Sakshi Gautam'' ^ 

Handloom street, H.No.-B-9-1076,Surgapuri,Kotkapura (Distt. Faridkot)-151204, India 
^Department of Physics, Panjab University, Chandigarh -160 014, India. 

In the present work, we make confrontation of our theoretical calculations using quantum 
molecular dynamics model with the experimental data for the reactions of ^°v4r +^^Sc, 
^^"^Au +^^^Au and ^"^^Xe +^^^Sn at different incident energies. In these reactions, we 
display the charge distribution and energy dependence of fragments multiplicity. Our 
results indicate good agreement with the experimental data for all the reactions . 



Email: sakshigautni@gniail.coin 



1 



1 Introduction 



Nuclear physics, in general and heavy-ion collisions, in particular, are of central interest 
due to several rare phenomena emerging at different incident energies. These include (1) 
collective flow and its disappearance [1], (2) breaking of colliding nuclei into pieces i.e. 
multi-fragmentation [2,3], (3) sub-threshold particle production [4] and (4) formation of 
hot and dense nuclear matter [5]. At high excitation energy, the colliding nuclei may break- 
up into several small and intermediate size fragments along with large number of nucleons 
that are also emitted. This is known as multi-fragmentation. The multi-fragmentation is 
the complex phenomena in nature which depends crucially on the incident energy of the 
projectile and on the geometry of the reaction [6]. At very low excitation energies, the 
excitation energy deposited in the system is too small to allow the break up of the nuclei 
into fragments. With the increase in the incident energy, the nuclei, after collision, can 
break up into dozens of fragments consisting of light, medium and heavy fragments. This 
phenomena was observed with the help of accelerators by Pcrfilov et al. and Lozkin et al. 
[7,8]. After 1982, Jakobsson et al. [9] observed a multiple emission of fragments of medium 
mass known as intermediate mass fragments (IMF's) in emulsion irradiation by the carbon 
beam of 250 MeV/nucleon at Berkeley BEVALAC. This scenario of multi-fragmentation 
was discussed in 1983 by Siemens [10]. Actually, the nuclear multi-fragmentation was 
discovered through cosmic rays accompanying the collisions of relativistic protons with 
targets and following the emission of slow fragments. These fragments were heavier than 
alpha particles but lighter than fission fragments. Also, the rise and fall of IMFs has 
also been established by Ogilvie et al. [11] through the collisions of Au particles with C, 
Al and Cu targets using ALADIN forward spectrometer at GSI, Darmstadt, with beam 
accelerated by SIS synchrotron. 

The development of accelerators led to huge experimental data on multifragmentation. 
On theoretical front also, a large number of models have been proposed to study multi- 
fragmentation. These include statistical models and dynamical models. In the present 
paper, we aim to show the validity of quantum molecular dynamics model by confronting 
the calculations with the experimental data over wide range of energies. 
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2 The Formalism 



2.1 Quantum Molecular Dynamics (QMD) model 

We describe the time evolution of a heavy-ion reaction within the framework of Quan- 
tum Molecular Dynamics (QMD) model [12-15] which is based on a molecular dynamics 
picture. Here each nucleon is represented by a coherent state of the form 
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where L is the Gaussian width of the particle. Thus, the wave function has two time 
dependent parameters Xa and Pa- The total n-body wave function is assumed to be a 
direct product of coherent states: 

= (l>a{Xl,Xa,Pa,t)(f)fs{x2,XfS,Pf3,t)...., (2) 

where antisymmetrization is neglected. One should, however, keep in the mind that the 
Pauli principle, which is very important at low incident energies, has been taken into 
account. The initial values of the parameters arc chosen in a way that the ensemble 
{At+Ap) nucleons give a proper density distribution as well as a proper momentum 
distribution of the projectile and target nuclei. The time evolution of the system is 
calculated using the generalized variational principle. We start out from the action 



r C[ct>A*]dT. (3) 
with the Lagrange functional 
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dt 



(4) 



where the total time derivative includes the derivatives with respect to the parameters. 
The time evolution is obtained by the requirement that the action is stationary under the 
allowed variation of the wave function 

5S^5 r C[(t),(t)*]dt^Q. (5) 

If the true solution of the Schrodinger equation is contained in the restricted set of 
wave function (j)^ {xi^XajPo) , this variation of the action gives the exact solution of the 



Schrodinger equation. If the parameter space is too restricted, we obtain that wave func- 
tion in the restricted parameter space which comes close to the solution of the Schrodinger 
equation. Performing the variation with the test wave function (2), we obtain for each 
parameter A an Euler-Lagrange equation; 

ddC dC 

For each coherent state and a Hamiltonian of the form, 

H = Ta + ^ ' Lagrangian and the Euler-Lagrange function can be 

easily calculated [15] 



^ 2Lm' 



m ^ 

Pa = -Vx„ (9) 

Thus, the variational approach has reduced the n-body Schrodinger equation to a set of 6n- 
different equations for the parameters which can be solved numerically. If one inspects the 
formalism carefully, one finds that the interaction potential which is actually the Briickner 
G-matrix can be divided into two parts: (i) a real part and (ii) an imaginary part. The 
real part of the potential acts like a potential whereas imaginary part is proportional to 
the cross section. 

In the present model, interaction potential comprises of the following terms: 

Va^ = m + Vil + Vcoul + VYuk, (10) 

Vioc is the Skyrme force whereas Vcoui and Vvuk define, respectively, the Coulomb and 
Yukawa potentials. Yukawa term separates surface which also play role in low energy pro- 
cess like fusion and cluster radioactivity [16,17]. The expectation value of these potentials 
is calculated as 

X(fr^d^rpd^p^d^p^, (11) 
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Vllc = y /a(Pa,ra,t)/^(P/3,r^,t)/7(P7,r^,t) 

x£pad^PfS(fp^. (12) 

where /a(Pa, i"a) ^) is the Wigner density which corresponds to the wave functions (eq. 2). 
If we deal with the local Skyrme force only, we get 



At+Ap 
-ySkyrme ^ 

a=l 



(13) 



Here A, B and C are the Skyrme parameters which are defined according to the ground state 
properties of a nucleus. Different values of C lead to different equations of state. A larger value 
of C (= 380 MeV) is often dubbed as stiff equation of state.The finite range Yukawa (Vyuk) and 
effective Coulomb potential (Vcoul) read as: 

Vy., = t3^^^4^-^f^, (14) 



The Yukawa interaction (with t3= -6.66 MeV and /j, = 1.5 fm) is essential for the surface 
effects. Also, the relativistic effect does not play role in low incident energy of present interest 
[18]. 

The phase space of the nucleons is stored at several time steps. The QMD model does not 
give any information about the fragments observed at the final stage of the reaction. In order 
to construct fragments from the present phase-space, one needs the clusterization algorithms. 
So, to construct fragments, we use Minimum Spanning Tree (MST) method. 
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Figure 1: The charge distribution (dn/dz) for the central coUision of ^°Ar +^^Sc at an 
incident energy of 85 MeV/nucleon. The experimental data has been taken by the Ref. 
[21]. 

3 Minimum spanning tree (MST) method 

The widely used clusterization algorithm is the Minimum Spanning Tree (MST) method [19]. In 

MST method, two nucleons are allowed to share the same fragment if their centroids are closer 
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than a distance rmin, 

\ri - < rmin- (16) 

where ri and rj are the spatial positions of both nucleons. The value of r^j„ can vary between 2-4 
fm. This method gives a big fragment at high density which splits into several light and medium 
mass fragments after several hundred fm/c. This procedure gives same fragment pattern for 
times later than 200 fm/c, but cannot be used for earlier times. 

We simulate thousand events for the reactions of ^^Ar +^^Sc (at 85 and 115 MeV/nucleon), 
^^'^Au +^^^Au ( at 150 MeV/nucleon) and ^^SXe+^^^Sn (at 25, 50 MeV/nucleon). The impact 
parameters have been guided by the corresponding experimental data. We use a soft equation of 
state along with standard energy-dependent nucleon-nucleon Cugnon cross-section. The details 
of the cross-section has been found in Ref. In fig. 1, we display the reaction of ^°Ar -|-^^Sc 
at central collision at an incident energy of 85 MeV/nucleon. Here, we display the charge 
distribution (open circles) for this reaction. The figure shows a linear decrease in the value of 
charge distribution (dn/dz) with the charge. This negative slope of charge distribution indicates 
a gradual transition from the spectator matter to the disassembly of the system. 

In the figure 2, we display the charge distribution for the reaction of "^^Ar +^^Sc at central 
collision at an incident energy of 115 MeV/nucleon. Here also, a decrease in the value of charge 
distribution is observed. Also, when the beam energy is increased from 85 MeV/nucleon to 115 
MeV/nucleon, the slope is still steeper indicates the total disassembly of the matter. When we 
compare our theoretical results with experimental data, we see that our theoretical calculations 
matches well with the experimental data at both the energies. This experimental data are taken 
with Michigan State University (MSU) FOPI array [20,21] at the National Superconducting 
Cyclotron Laboratory (NSCL) using beams from K1200 Cyclotron. The FOPI array consists of 
a main ball of 170 phoswich counters covering angles from 23 degrees to 157 degrees. 
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Figure 2: Same as figure 1, but at an incident energy of 115 MeV/nucleon. The experi- 
mental data has been taken by the Ref. [21]. 

In fig. 3, we display the results for the central colUsions of ^^'^Au +^^^Au colUsions. This 
reaction was measured at the incident energy of E=150 MeV/nucleon. Also, the condition of 
rapidity cut is applied so that there may be the covering of only those particles which are emitted 
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Figure 3: The multiplicity distribution (dn/dz) of fragments measured for ^"^"^ Au +^^^Au 
reaction at 150 MeV/nucleon having central collisions. The experimental data has been 
taken by the Ref. [23,24]. 

from the most violent reactions. This cut also removes the spectator particles i.e. only those 
particles are taken which took part in the reaction. 
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Figure 4: The charge distribution for central collisions of ^^^Xe+^^^Sn at bombarding 
energy of 25 MeV/nucleon. The experimental data has been taken by the Ref. [26,27]. 

In the figure, we again see a decrease in the value of charge distribution (dn/dz) with the 
increase in charge. The reason is same as mentioned earlier in the description of figure 1. In 
the figure, the open circles show our theoretical calculations. When we compare our theoretical 
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Figure 5: Same as Figure 4, but at incident energy of E= 50 MeV/nucleon. The experi- 
mental data has been taken by the Ref. [26,27]. 

calculations with the experimental data, it matches well with each other. These experimental 
data are taken with the phase 1 set up of FOPI facility. The Au ion beam was delivered by the 
rapid-cycling synchrotron SIS, Darmstadt [22-24]. 
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Figure 6: The IMF multiplicity versus beam energy is shown for central collision in the 
reaction of ^"^^Xe +^^^Sn. The experimental data has been taken by the Ref. [28]. 

In fig. 4, we display the central collisions (b/b^ax=0.0) of ^"^^Xe +^^^Sn at bombarding 
energy of E=25 MeV/nucleon. Open circles shows our theoretical calculations done by the MST 
approach. In compare to previous results, we see now two separate zones. First, a negative 
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slope for light fragments and 'u' shape for heavy fragments. This indicates that at low energy, 
reaction produces excited compound system, thus cool down with the emission of free nucleons. 
Nearly, no intermediate mass fragment is seen in the reaction. 

In figure 5, we again display the central collisions (b/bmaa;=0.0) of ^^^Xe +-'^^^Sn at bom- 
barding energy of E= 50 MeV/nucleon. We again see a similar steepening of charge distribution 
with the charge. When we compare our theoretical calculations with the experimental data 
taken by the INDRA at the GANIL and GSI accelerator [25-27], they are in good agreement 
with each other. 

In figure 6, we display the reaction of ^^^Xe -|-^^^Sn reaction for central collisions. We take 
only those particles with Z^3. In the figure, we have seen the variation of number of intermediate 
mass fragments with the increase in the beam energy. At low energy i.e. 25 MeV/nucleon, 
number of IMFs produced is less because we have mainly the heavy fragment and free nucleons. 
Also, there is more time for the reaction to take place at low energy as compared to high energy. 
But with the increase in energy, there is increase in IMFs as shown in the figure. When our 
theoretical calculations are matched with the experimental data taken during the fragmentation 
studies performed with FOPI detector [28], they are in good agreement. 

Summary 

We studied the different symmetric and asymmetric reactions mainly for central collisions with 
different range of energies. We observed the steepening of charge distribution with the increase 
in charge. Also, we observed that number of IMFs increases with the increase in beam energy. 
Our calculations showed good agreement with the data at all the energies. 
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